Improved Unitarized Heavy Baryon Chiral Perturbation Theory for tyN Scattering to 

fourth order 

A. Gomez Nicola 1 ^ J. Nieves 2 Q J.R. Pelaez 1 !]] and E. Ruiz Arriolafl 

1 Departamento de Fisica Tedrica. Universidad Complutense. 28040 Madrid, Spain, and 
2 Departamento de Fisica Moderna, Universidad de Granada, E-18071 Granada, Spain 

We extend our previous analysis of the unitarized pion-nucleon scattering amplitude including up 
to fourth order terms in Heavy Baryon Chiral Perturbation Theory. We pay special attention to the 
stability of the generated A(1232) resonance, the convergence problems and the power counting of 
the chiral parameters. 
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I. INTRODUCTION 



Unitarization methods have been widely and successfully employed in the recent past to enlarge the applicability 
region of Chiral Perturbation Theory (ChPT) expansions, both in the meson-meson sector as well as in the meson 
baryon sector and to describe the lightest resonances without including them explicitly as degrees of freedom. Two 
important constraints are required: exact unitarity and compliance with Chiral Perturbation Theory at a given order 
' of the expansion. In practice this approach provides a remarkable description of data in the scattering region. In 
the case of ttN scattering in the elastic region, the subject of this paper, a thorough partial wave analysis exists Q 
(see also the recent update 0)- For such a system, pions and nucleons are treated as explicit degrees of freedom and 
a consistent counting becomes possible if nucleons are treated as heavy particles but in a covariant framework 
yielding the so called Heavy-Baryon Chiral Perturbation Theory (HBChPT) QUE El- In this counting, the expansion 
for the scattering amplitude is done as a series of e N / (F 21 M N+1 ~ 21 ) terms , with I = 1, . . . , [(N + l)/2], M the 
baryon mass, and F the pion decay constant. The quantity e is a generic parameter with dimensions of energy 
■ constructed in terms of the pseudoscalar momenta and the velocity (v 2 = 1 ) and off-shellness k of the baryons 

defined through the equation pb =M v + k, with ps the baryon four momentum and M the baryon mass at leading 
order in the expansion. After the relevant effective Lagrangian was written down 0, and the issue of wave function 
renormalization was studied Q standard HBChPT calculations to second third |ToL ITlj and fourth ^| order 
have become available. The unitarization of these amplitudes of irN scattering in the elastic region has followed 
closely these developments, particularly the third order calculation [lollllj . This is the lowest order which generates a 
perturbative unitarity correction of the amplitude. The unitarization was carried out either using the standard Inverse 
Amplitude Method [u| (IAM) or its improved version pj| 1 . By successful we mean the possibility of describing the 
data in the resonance region with parameters of natural size. The purpose of the present paper is to extend the study 
initiated in Ref. [l4| and to analyze specifically the qualitative and quantitative new effects generated by the fourth 
order contribution calculated in Ref. |l2( in our unitarization scheme. 

Let us then specify the scope and motivations of our work: First, our scheme is based on two fundamental ideas: 
demanding exact unitarity and considering the F~ 2 HBChPT expansion independent from and converging faster than 
the M _1 one. In 0] we showed that this allows to generate the A(1232) as well as to fit the remaining S and P wave 
channels with natural values for the low-energy constants (LEC) unlike for instance the IAM 13]. Our method was 
implemented in pdj with the first contribution of order F~ 4 only, coming from the third order amplitude. Including 
the fourth order will allow us to check the convergence of our method by considering, for instance, the 0(F^ 4 M~ 1 ), 
to be included in the third order F~ 4 term. 

Second, there is an interesting issue that we did not account for in 01 which has to do with the separation of the 
dimensionful third and fourth order LEC into two pieces contributing to the orders F~ 2 and F~ 4 . As we will see 
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below, taking into account this effect may change considerably our description of the partial waves. The reason why 
we did not consider it in [l4| is that we used the amplitudes in |10| . which provide an specific separation that turns 
out to be very natural, as we will see below 2 . 

Third, comparing the perturbative results to order three [Til ] and four |T^. one observes that in order to achieve a 
reasonable convergence, the fourth order constants become of unnatural size and, furthermore, their particular values 
are often incompatible from one fit to another. This is a signal of the bad convergence of the HBChPT series and 
could influence also the convergence of our unitarized formula. 

Fourth, unitarization methods are rarely applied beyond the leading order in the imaginary part of the amplitudes 
[lii Il7j. The study of the fourth order of ttN system within HBChPT provides an opportunity to learn about the 
unitarization approach beyond this lowest order. 

For comparison purposes with previous works ITol ITU IT2I Il3l Il4| we will take the partial wave analysis performed 
in Ref. Q]. The recent update f| does not bring significant changes to our discussion. 



II. THE UNITARIZED AMPLITUDE 



In order to have a neat separate expansion of the partial waves in powers of M~ x and F~ 2 , we need to re-expand 
the amplitudes in [l^ . as it was already done to third order in ^3] with those in Then, following the notation 
in [T^j. we have, to fourth order, for any partial wave 
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with m the pion mass, M the nucleon mass, F the pion decay constant and to the pion CM energy. The partial wave 
unitarity condition 



Im/^ 1 = -q, 

where q is the CM momentum, implies that perturbatively one has 3 
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Following the same ideas as in Il4| . we will consider the unitarized amplitude to fourth order 
1 F 2 1 



£(3,3) + l|t( 3 < 4 ) 



/lunitarizcd " m *(!.!) + S f(l,2) + ( ™ )2 <(1,3) + ( m )3 f (l,4) [t( l,l)]2 + 2 m t (l,l) t (l,2) 

which, using @ yields immediately (J2J). Let us recall that our Improved IAM formula at third order read ^ 



F 2 



f\ 



Unitarized 



4(1,1)+ ™ t (l,2)+ (S)2 t(l,3) 
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A similar situation has appeared already in the NNLO uni tar y analysis of tttt scattering Il7l . 

We have checked analytically that the amplitudes in are perturbatively unitary if the following misprints are corrected: 
^7 — M should read —^—^-(-uj) in their eq.(3.16) and 6u) 4 {-4Ml + 4u) 2 + t) should read -6lu 4 {-AMI + 4ui 2 + t) in their 

eq.(3.18) . In fact, with these two signs corrected, we reproduce the threshold parameter expressions given in their eqs.(A.l)-(A.8), except 
for the 7r 2 in the denominator of the fourth term in the r.h.s. of their (A. 8) which should read 7r 3 and the +<7a<?18 M 2 / (4ir F 2 (M n +m)) 
in eq.(A.3), that should have the opposite sign. 



3 



which can be now reobtained from eq.Q by removing the tl 1 ' 4 ) and i' 3 ' 4 ) terms and, consistently with unitarity, 
removing also the 2(m/M)t^ 1 ' 1 't^ 1 ' 2 ' factor in the second denominator. Hence, as we have stressed in the introduction, 
the knowledge of t^ 1,4 ^ and ^ 3 ' 4 ^ allows us to test our power counting by including one more term both in the 0(F~ 2 ) 
and 0(F~ 4 ) contributions. 

III. THE THIRD ORDER AND THE LEC POWER COUNTING 

In the literature there are two 0(q 3 ) calculation s TToL ITU . using different choices of counterterms and renormalization 
schemes, but only one at 0{q i ) ljj] following the |ll| scheme. The translation between them does not simply amount 
to a change of notation, but involves some 1/M corrections. Since our results at third order |l4j | were constructed 
directly from yLfJ, we have to check to what extend our previous 0(q 3 ) results are reproduced when using the 
amplitudes and notation of Refs. [Tll[T^| . In so doing two remarks are in order: 

First, already at third order, the re-expanded amplitudes of ^(J an d ^3 differ slightly due both to a different 
choice of the reference frame and of the nucleon wave function renormalization (see comments in In practice, 

this just means that there are slight numerical differences between the perturbative results, which eventually could 
be absorbed in the numerical values of the LEC of the HBChPT Lagrangian. 

A second point becomes more relevant for our purposes: if we just take the third order amplitudes in J Hi. we 
re-expand them separating the different contributions and use our third order unitarized formula (eq.(13) in |l4|). we 
find a much worse result than in [l4|, particularly in the P33 channel where the A(1232) resonance should appear. 
This is shown in Figure^. We remark that we are using a set of parameters compatible with those in |14| . where 
the description of the resonance was excellent within the errors even without fitting. 

The origin of this apparent discrepancy is that we have not taken into account that in our power counting scheme, 
the LEC themselves may have contributions of different orders. In fact, all the difference with 0] is that we have 
chosen now a different parametrization of LEC, although the numerical values are compatible: ™ we followed 
Vn. [LOB where the five 0(q 3 ) LEC appearing in the amplitude are called bi + b 2 ,b 3 ,b 6 , 6 16 — 6 15 , 6 19 . Here, we follow 
[ill If^ l where the relevant 0(q 3 ) constants are d\ + di , J3 , J5 , J14 — Ji5,Ji8. Comparing the Lagrangian given in 
equations (2.45)-(2.47) of [HI with that in [7j one observes that the bi are related to the di for i = 1, 2, 3 and to the 
di-i for i = 6, 15, 16, 19 typically as: 

M 2 di ~ constant + b t M 2 /(Wir 2 F 2 ) (6) 

with a constant that is 0(1 ) in the F~ 2 counting. This comes from the fact that in some finite terms coming from 
renormalization have been absorbed in the bi. 

Now, following our power counting arguments, if we replace in the amplitudes of refs.jlll Il2j| the di using eq.jfjj), 
there are pieces in t' l,3 > shifted to £( 3,3 ) (remember that all the dependence with the di is in t^ 1 ' 3 )). This changes 
the functional dependence of t^- 3 ' 3 \ including, for instance, higher order polynomial contributions that otherwise were 
not present. Although the perturbative amplitude remains the same, the unitarized one changes since i' 1 ' 3 ) and £( 3 ' 3 ' 
are treated on a different footing. With this procedure we obtain the unitarized results shown in Figure ^3. The 
improvement is clear for the P33 wave and the results are similar to those in [l4|. The corresponding values for the 
mass and width of the A(1232) extracted from the phase shifts are given in the second column of Table [I] This 
highlights the importance of taking into account the counting of the LEC. 

The above separation is, of course, arbitrary, since nothing prevents us from normalizing the di, which are quantities 
of dimension E~ 2 as (AnF 2 )di instead of M 2 di, assuming that both are quantities of natural size. Thus, the most 
general way to proceed would be to consider as free parameters the coefficients of the 0(1) and (D(F~ 2 ) terms in M 2 di. 
In such a way we would duplicate the number of 0(<z 3 ) LEC, rendering the approach unnecessarily complicated, since 
we already know that it is enough to consider the separation given by the bi parametrization |2J. We will thus use 
only that separation in our calculations, but, after performing the fit, we will give the results in the di set for easier 
comparison with the literature. 

The results of our 0(g 3 ) fit is shown in Figure|2] The fit parameters and their errors are given in the fourth column 
of Tablc lLTTl whereas the results for the A mass and width are given in TableHl All them are in agreement with what we 
found in [14(- The description of data in all channels is very good and our 0(<7 3 ) LEC are all of natural size although, 
as it also happened in [l4|. they differ somewhat from those obtained from HBChPT (second column in Table ITTT1 . 
Note that systematic errors are not given in Table III, although they are dominant as it can be seen from Table II). 
This is probably another consequence of the poor convergence of the HBChPT series. Let us nevertheless recall that, 
in 0] it was shown that one can perform 0(q 3 ) fits where the 0{q 2 ) parameters Cj are fixed to the predictions of 
Resonance Saturation [!j and the results are still in excellent agreement with data. 
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FIG. 1: 0{q i ) Unitarized phase shifts as a function of the total CM energy y/s: As explained in section III: a) In the two left columns 
we use the <i;, whereas in b) in the two right columns we use the bj s et. Experimental data are from Q. The areas between dotted lines 
correspond to the propagated errors of the parameters of fit 1 in llll in both cases. 
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IV. FOURTH ORDER RESULTS 

Let us consider now our fourth order unitarized amplitude (0} with the HBChPT amplitudes of I n principle, 
the 0(q 4 ) amplitude depends on nine different combinations of 0(q 4 ) constants e,, in addition to the four 0(q 2 ) Ci 
and the five 0(q 3 ) dj. These nine combinations are displayed in the first column of Table tTTTl However, as noted in 
the last four combinations actually amount to renormalize the c,, giving rise to new c, as given in eq.(3.23) of 
|L2|. Strictly speaking, the 0(q 4 ) amplitude depends on 18 parameters, since the still appear in the pure 0(q 4 ) 
terms. However, replacing those Ci by Cj introduces higher order corrections, so that the number of free parameters 
up to 0(q 4 ) is really 14. At this point, as commented in ^1 one can follow two different strategies. The first one is to 
consider as free parameters Ci,di and the with i = 14 — 18. This is the parameter set listed in Table ITT1 Although 
this is the more natural set, it has the inconvenience that one cannot disentangle which part of the Ci comes from 
0(q 4 ) renormalization, since those corrections are relatively large (another signal of the HBChPT bad convergence). 
As a consequence, it becomes more difficult to compare with previously published values for the c$. The alternative 
(strategy 2) is to fix the Ci values, which in turn are the ones less subjected to uncertainties and then use the di 
and the nine combinations of as free parameters. This second strategy is useful for instance to fix the Ci to the 
predictions of Resonance Saturation || as we also did in ^4) . 

In addition, we have to face again the problem of the LEC counting, according to the discussion in the previous 
section. Thus, besides the "reordering" of terms coming from the di, we now have to consider also that coming from 
the separation of M 3 e~i into a constant (a contribution to v- 1 '*') plus an O (M 2 /(4ttF) 2 ) term (which contributes to 
i( 3 - 4 )). Recall that in this case we do not have any "natural" way to perform that separation, as in the 0(q 3 ) case. 

A. The unitarized partial waves to 0(q 4 ) 

First, as we did to 0(q 3 ), we will show the predictions of our formula without fitting, performing a Monte Carlo 
sampling of the perturbative LEC, assuming that they are uncorrelated. Following the first strategy, we have used the 
LEC given in in particular those given in their "Fit 3" that we reproduce in Table |Hj We also list the "Fit 1 and 
2" parameter sets to illustrate that, as pointed out in |l2j, the systematic errors are much larger than the statistical 
ones, that we are quoting in the Table. For that reason we will take bigger errors, since the errors listed in are 
clearly underestimated. In view of the uncertainties in [l2T | we have assigned an error of 1.0 to ei4, eis, en, e^, of 0.5 
to ci, C2, C3 (the Ci have bigger uncertainties than the Ci due to their contribution) and of 0.25 to the remaining 
LEC. In Figure^ we show the 0(q 4 ) prediction "redefining" the di as before but without doing so for the e^, while 
in Figure^ we have also redefined the e~i for convenience as M 3 e.i = 1 + fiM 2 /(4ttF) 2 . Throughout this paper, and 
for practical purposes, we will consider only these two situations. 

In view of Figure |21 there are several comments in order: First, consider the P33 channel, where our approach is 
meant to be more accurate. Here our 0(q 4 ) result confirms the 0(q 3 ) one and even improves it slightly. Observe for 
instance the results for the A parameters given in the fourth column of Table |U corresponding to Figure 03. This is 
one of the main conclusions of this work, namely that the 0(q 4 ) calculation confirms that our unitarization method 
generates dynamically the A(1232) resonance. The improvement of the P33 channel description is also a common 
feature with [l^. Note also that this conclusion does not change by using the or the fi formulas as long as we 
redefine the 0{q 3 ) LEC as before. 

As for the other channels, we see no actual improvement when comparing with the unfitted 0(q 3 ) in Figure \T]p. 
On the contrary, we get worse results for most of them, especially the two S channels and the Pn one. Here, we see 
significant differences between using the fi prescription or not. In fact, without fitting, our choice for the fi does not 
seem to give better results than using the gj directly (Figure^) or, equivalently, neglecting the (D(F~ 2 ) contribution 
in M 3 Ei. 

The hope is that we can perform 0{q ) fits which improve these five channels without spoiling the P33 one and 
with a reasonable size for the LEC. However, we must bear in mind that, as commented before, the low-energy 0(q 4 ) 
fits performed in ^3 already show that one gets only slightly better descriptions and bigger uncertainties for the LEC 
than the 0(q 3 ). 

B. 0(q 4 ) fits 

In Figure 0Ji we show the result of our best fit with the fit errors propagated. Here we have followed the first 
strategy and we have used the fi defined in the previous section. The LEC and their errors are given in the last 
column in Table UTI The main observation is that we reproduce the data with constants of natural size. The constant 
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FIG. 3: 0(q 4 ) Unitarized phase shifts as a function of the total CM energy y/s. a) The two left columns correspond to reordering the di 
only, whereas b) in the two right columns we have reordered both d~i and e^. Experimental data are from 0. The areas between dotted 
lines correspond to the propagated errors of the parameters of fit 3 in |12| with bigger errors (see Sect. IV. a) 





0(q 3 ) unfitted 


fit 0(q 3 ) 


0(q 4 ) unfitted 


fit 0(q 4 ) 


PDG 


M A (MeV) 


1221 t» 


1229 t\t 


123819° 


12321^ 


1230 - 1234 


T A (MeV) 


111 2+ 16 - 9 


108.4 ±™i 


125.2±^ 


107.3 ±£.1 


115-125 



TABLE I: A(1232) resonance parameters in the different cases considered in this paper, 
obtained from the condition Sl 3 \ s=M 2 = ty/2 and 1/Ta = MA(d8l 3 /ds)\ s=M 2 . 



The resonance mass and width are 



eis turns out to be highly correlated numerically with e~i7 so that we have chosen to fix one of them to the perturbative 
value. For the A parameters we get the results in the fifth column in Table [I] which is still fully compatible with the 
experimental result. Note that, as expected from our previous comments, the uncertainties in the fit parameters are 
now bigger than in the 0(q 3 ) fit. However, the quality of the fit is comparable if not better: we get a y 2 Jd.o. f ~ 0.3 
for the 0(q 3 ) fit in FigureEJand y^/d.o.f ~ 0.17 for that in Figure^. As is customarily done for the x 2 

calculation we havee added some error to the data; in particular, we have chosen to add a 3% relative error plus one 
degree systematic error. Therefore, our method shows clear signs of convergence when we perform unconstrained fits, 
although the uncertainties in the LEC remind us of the bad convergence of the HBChPT series. Note also that the 
bigger uncertainties are in the S channels, as we have commented before. 

We also show, in Figure the result of a fit using the second strategy and fixing the Ci as the central values of 
the predictions of Resonance Saturation As it also happened with the 0(q 3 ) in [14(, the fit result is slightly worse 
when the Ci are not free parameters. Here we obtain a better fit when using directly the e, and not the fi as free 0{q A ) 
parameters. Nevertheless, some of the constants become of unnatural size and also the errors are bigger than for the 
unconstrained fit. For the fit in Figure QJd we get a x 2 /d-0-f ~ 0.58 and the LEC listed in the fifth column in Table 
IIII The correlations between the different LEC also become more important. Here, in addition to en and eis, there 
are also strong correlations among en, eis, ei6, £22 — 4e"38, §20 + 635 — gAdi6/(8M), 2e~2i — ew and 2eig — £22 — e~36 
which allow to fix one of them. It should be commented that one could find fits with more natural values and a higher 
X 2 /d.o.f but we have preferred to show the best fit, emphasizing the convergence problems. 
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Fettes-Meifiner [12] 


Our IAM 




Fit 1 


Fit 2 


Fit 3 


0(q 4 ) fit 


Cl 


-2.54 ±0.03 


-0.27 ±001 


-3.31± 0.03 


- 1.43 ± 0.10 


C2 


0.60 ± 0.04 


3.29 ± 0.03 


0.13 ± 0.03 


- 0.33 ± 0.10 


C3 


-8.86 ± 0.06 


-1.44 ±0.03 


-10.37 ± 0.05 


- 2.62 ± 0.17 


C4 


2.80 ±0.13 


3.53 ± 0.08 


2.86±0.10 


0.64±0.10 


di + d,2 


5.68 ± 0.09 


4.45 ± 0.05 


5.59± 0.06 


1.11 ± 1.02 


d 3 


-4.82 ± 0.09 


-2.96 ±0.05 


-4.91±0.07 


-0.59±1.06 


d 5 


-0.09 ± 0.06 


-0.95 ± 0.03 


-0.15±0.05 


-0.19 ± 0.40 


dl4 — dl5 


-10.49 ± 0.18 


-7.02 ±0.11 


-11.14 ± 0.11 


2.49 ±1.93 


dia 


-1.53 ±0.17 


-0.97 ±0.11 


-0.85±0.06 


- 17.49± 1.31 


E14 


6.39 ± 0.27 


-4.68 ± 0.14 


7.83±0.23 


1.58 ± 0.53 


§15 


4.65 ± 0.31 


-18.41 ±0.15 


9.72±0.25 


- 1.41 ± 1.13 


§16 


7.05 ± 0.30 


7.79 ±0.15 


6.42± 0.25 


3.50 ± 1.30 


§17 


4.88 ± 0.98 


-17.79 ±0.53 


5.47± 0.64 


6.56 ± 1.92 


§18 


-9.15 ±0.98 


19.66 ±0.53 


-0.17±0.64 


-0.17 (fixed) 



TABLE II: HBChPT low energy constants with strategy 1. The Cj, di and are given in GeV 1 ,GeV 2 and GeV 3 respectively. We list 
the three parameter sets provided in \\ A to illustrate the large systematic uncertainties already existing in the perturbative determinations. 



V. CONCLUSIONS 



In this work we have used a unitarized irN scattering amplitude including up to 0(q 4 ) terms in the standard Heavy 
Baryon Chiral Perturbation Theory expansion. This has allowed to test our method of considering the F~ 2 expansion 
resumming the M contributions. The description of the P33 channel and the A resonance, which is dynamically 
generated, are excellent within the experimental errors. The inclusion of the new 0(q 4 ) terms does not change much 
this picture, which is a consistency check of our formalism. 



8 





Fettes et al. fill 


Fettes-Meifiner [12] 
0(q 4 ) (Strategy 2) 


0(q 3 ) fit 


0(q*) fit 
RS a of [9J 


Cl 


-1.53± 0.18 


-1.47 (input) 


-0.43 ± 0.04 


-0.9 (input) 


C2 


3.22 ± 0.25 


3.26 (input) 


1.28±0.03 


3.9 (input) 


C3 


-6.20± 0.09 


-6.14 (input) 


-3.10 ± 0.05 


— 5.3 (input) 


C4 


3.51±0.04 


3.50 (input) 


1.51±0.04 


3.7 (input) 


di ± d 2 


2.68 ±0.15 


4.90 ± 0.05 


2.66 ±0.20 


10.36 ±0.53 


da 


-3.11±0.79 


-4.19±0.07 


-0.32 ±0.2 


- 4.07 ±0.26 




0.43±0.49 


-0.16 ±0.05 


-1.66 ±0.10 


- 3.23 ± 0.31 


dl4 — di5 


-5.74 ±0.29 


-9.31 ±0.10 


-5.34 ±0.40 


-1.17 ± 0.68 




-0.83±0.06 


-0.84±0.06 


-2.60±0.20 


- 53.29 ± 3.37 


ei4 




4.19 ±0.23 




2.24 ± 0.94 


eis 




4.54 ± 0.25 




- 2.17 ± 2.14 


ei6 




2.74 ±0.24 




5.12 ± 1.26 


ei7 




7.20 ± 0.64 




- 0.95 ± 1.86 


eis 




- 3.36±0.64 




- 3.36 (fixed) 


e22 — 4e38 




27.72±0.74 




27.72 (fixed) 


e20 + e 3 5 - gAdie/(8M) 




-17.35±0.36 




- 70.11 ± 3.07 


2ei9 — e 2 2 — e~m 




-25.12 ± 0.69 




97.14 ± 10.00 


2e~2i — e-37 




-5.00±1.43 




17.64 ±10.7 



TABLE III: HBChPT low energy constants with strategy 2. The Ci, di and ei are given in GeV ^GeV 2 and GeV 3 respectively. 



In order to describe accurately the data in the other three P channels and the two S ones, one needs to ht the LEC. 
We have showed that it is possible to fit the six channels simultaneously with natural values for the LEC, although 
with considerably larger uncertainties than in the 0(q 3 ). This is a consequence of the poor HBChPT convergence 
which shows up already at the pure perturbative level at lower energies. In fact, when one tries to perform 0(q 4 ) fits 
constraining the lowest order constants to the Resonance Saturation hypothesis, some of the LEC become of unnatural 
size and their errors increase considerably. These convergence problems are especially important in the two S waves. 

We have also discussed the issue of the LEC power counting, which is relevant in our expansion scheme. The 
importance of this effect has been especially highlighted in the P33 channel, where a correct splitting of the 0(q 3 ) 
constants is crucial. To 0(q 4 ) the influence of such LEC reordering is smaller as far as the resonance is concerned, 
although it may improve the convergence in the other channels. 

In summary, our unitarization method is robust and is almost not affected by the HBChPT convergence problems as 
far as the generation of the A(1232) resonance is concerned. However, the predictions of the unitarized amplitude to 
fourth order for other channels show similar problems of convergence as the perturbative one, even though one can still 
find excellent descriptions of data with natural values for the low-energy constants. It seems a natural continuation 
of this work to implement our unitarization methods within the context of the Lorentz invariant formalism proposed 
in [lj| (see also Ref. [lj|) which has better convergence properties. 
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